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By means of a simple rescaling, modifications of hadron masses and widths are 
incorporated into the thermal analysis of particle ratios in ultra-relativistic heavy- 
ion collisions. We find that moderate, up to 20%, changes of hadron masses do not 
(Si ■ spoil the quality of the fits, which remain as good as those obtained without mod- 

ifications. The fits with the modified masses yield, however, modified values of the 
optimal temperature and baryon chemical potential. In particular, with decreas- 
ing masses of all hadrons (except for pseudo-Goldstone bosons) the fitted values of 
the temperature and the baryon chemical potential are lowered, with the change 
approximately proportional to the scaling of masses. In addition, we find that the 
broadening of the hadron widths by less than a factor of two practically does not 
affect the fits. 
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Recent theoretical studies |],0||[|,|5|] show that the hadronic yields and ratios 
in ultra-relativistic heavy-ion collisions can be described well in the frame- 
work of simple thermal models. The thermodynamic parameters obtained in 
this approach define the so-called chemical freeze- out point, i.e., a stage in the 
evolution of the hadronic system when the "chemical composition" is estab- 
lished. Thermal-model fits to the SPS data show that the temperature at the 
chemical freeze-out, T chem , as well as the baryon density, are large. One typi- 
cally obtains T c h em ~ 170 MeV, which is close to the expected critical value for 
the deconfmement /hadronization phase transition. In this situation one may 
expect that hadron properties at the chemical freeze-out are strongly modified 
by the presence of the hadronic environment. Indeed, such modifications are 
predicted by model calculations Pffl|5|P, |Tu||lT] , |l3?[ , which inter alia helps to 
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explain the low-mass dilepton enhancement observed in the CERES [13] and 



HELIOS |14fl experiments. In this paper we incorporate possible modifications 
of hadron masses and widths into thermal analysis of particle ratios. We 
generalize the results of Refs. |15|J16[| where the problem was studied without 



refitting thermodynamical parameters. 

In the first part we include only the mass modifications and calculate the 
particle densities from the ideal-gas expression^ 
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p 2 dp 



27T 2 J exp [(E* fI^ hem Bi ^chem^i AherrJi) /T c hem] ± 1 

where <?j is the spin degeneracy factor of the ith hadron, Bi,Si,Ii are the 
baryon number, strangeness, and the third component of isospin, and E* = 



p 2 + (m*) is the energy. The quantities //f ftem , \i s chem and ^ chem are the 
chemical potentials enforcing the appropriate conservation laws. In standard 
thermal-model fits Eq. ([I]) is used with the vacuum masses, m* = m^. The in- 
medium masses, m*, may depend on temperature and density in a complicated 
way. In order to explore possible different behavior of in-medium masses and, 
at the same time, keep simplicity, we do our calculations with the meson and 
baryon masses rescaled by the two universal parameters, xm and xb, namely 

m* M = x M m M , rn* B = x B m B . (2) 

An exception from this rule are the masses of pseudo-Goldstone bosons (tt, K, rf) 
which we keep constant. 

Equation ([[]) is used to calculate the "primordial" density of stable hadrons 
and resonances at the chemical freeze-out. The final (observed) multiplicities 
receive contributions from the primordial stable hadrons, and from the sec- 
ondary hadrons produced by decays of resonances after the freeze-out. We 
include all light-flavor hadrons listed in the newest review of particle physics 
I9|, with a few exceptions for the cases where the properties of a listed par- 



ticle are ambiguous or not known sufficiently well. Ref. |nj] is also used to 
determine the branching ratios. We neglect the finite-size and excluded vol- 
ume corrections: the former are negligible, whereas the latter do not affect the 
particle ratios.0 



1 The use of Eq. (|l]) is valid when the in-medium hadrons can be regarded as good 
quasi-particles. A thermodynamically consistent approach has been constructed so 
far only for the lowest multiplets of hadrons |p^JI^| . At SPS energies, however, it 
is crucial to include all hadrons with masses up to (at least) 1.8 GeV. For such a 
complicated system, a thermodynamically consistent approach is not at hand at the 
moment. 

2 The last property is due to the fact that the overwhelming majority of hadrons 
is heavy and may be treated as classical particles. In this case we may replace the 
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Fig. 1. Dependence of x 2 > and the fitted values of the temperature and the baryon 
density on the scale parameter x. The plot is for the Pb + Pb collisions at SPS 
energies, with 19 particle ratios included in the analysis. Solid lines: all hadron 
masses (except for Goldstone bosons) are scaled with xm = xb = x. Dashed lines: 
only meson masses are scaled, xu = x,xb = 1- Dotted lines: only baryon masses 
are scaled, xb = x,xm = 1- The nuclear saturation density po = 0.17fm~ 3 . 

For given values of % and xb we fit the temperature, T c h e m, and the bary- 
onic chemical potential, /i^ em , by the minimization the expression x 2 — 



^ / r>exp Dtherm\' 



where R^ xp is the kth measured ratio, is the 
corresponding error, and _R* herm is the same ratio as determined from the 
thermal model. The strangeness chemical potential, \i s chem -, and the isospin 
chemical potential, /4/iem' are determined from the requirements that: i) the 



Fermi-Dirac (Bose-Einstein) distribution function in (|l|) by the Boltzmann distribu- 
tion function. Thus, for equal excluded volumes of all hadrons the excluded volume 
corrections factorize and cancel out in the ratios. We have checked that the use of 
the classical statistics changes our results by a few percent only. 
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initial strangeness of the system is zero, and ii) the ratio of the baryon number 
to the electric charge is the same as in the colliding nuclei. 



In Fig. 1 we plot our results obtained for the experimental ratios for Pb + Pb 
collisions at SPS, as compiled in Ref. ||. In the case xm = = 1 we find: 
T chem = 169 MeV, ^ hem = 250 MeV, /if hem = 65 MeV, and fx J chem = -9 MeV. 
These values are in good agreement with those found in Ref. 0: T c h em = 168 
MeV, /4 em = 266 MeV, /4 em = 71 MeV, and ^ chem = -5 MeV. In Fig. 1 a) 
we give our values of \ 2 . One can observe that a small decrease of the meson 
and baryon masses, % = xb ~ 0.9, leads to a slightly better fit with the 
corresponding smaller values of the temperature and the baryon density, as 
shown in Figs. 1 b) and 1 c). To check the significance of this result we have 
also analyzed Si + Au collisions at AGS, and S + Au collisions at SPS. In 
these two cases x 2 has a flat minimum at xm ~ %b ^ 1- We thus conclude 
that moderate modifications of hadron masses, say by 20%, do not spoil the 
quality of thermal fits. However, with the modified masses the thermodynamic 
parameters characterizing the fits do change. For example, if we rescale both 
meson and baryon masses (except for Goldstone bosons) in the same way, 
x = xm = xb, the temperature and the chemical potentials are to a very 
good approximation also rescaled by x. This follows from the fact that we 
study a system of equations which is invariant under rescaling of all quantities 
with the dimension of energy. If we allowed also for the changes of the masses 
of the Goldstone bosons, the thermodynamic parameters would scale exactly 
as T chem (x) = x T chem (x = 0) and ^ c hem{x) = x fi c hem{x = 0). In this case \ 2 
remains constant, independently of x. For fixed values of the Goldstone-boson 
masses, the scale invariance is broken, and x 2 varies weakly with x, as shown 
in Fig. 1 a). 

To account for finite in-medium widths, T*, of the resonances we generalize 
Eq. (D to the formula |(|gl],El,il 
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rii = J dM 2 J dp 
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27T 2 exp 



\fM 2 + p 2 Hchem-Bi ^chem^i ^chem^i) /^i 



- chem 
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where N is the normalization of the relativistic Breit-Wigner function, N = 
\ + - arctan [((m*) 2 — M%) / (m*T*)\ I. The integral over M 2 is taken to 
start at the threshold Mq corresponding to the dominant decay channel. In 
the limit r* — > Eq. @ obviously reduces to formula (p. 

In order to analyze the effect of broadening of hadron widths we introduce the 
parameter y in such a way that 

I7 = yIV (4) 
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Fig. 2. Dependence of x 2 , and the optimal values of the temperature and the baryon 
density on the scale parameter for the widths, y. 



Here Tj are the vacuum widths, hence the case y — 1 corresponds to the 
physical widths as measured in the vacuum, and the case y = represents the 
situation when the widths are neglected (our previous analysis based on Eq. 
([!])). In Fig. 2 we show the results of our fitting procedure. We observe that 
the inclusion of the vacuum widths does not change the value of x 2 ; and the 
values of T c h em and p® hem - An increase of the widths by a factor of 2 has also 
little effect. Only for larger modifications of the widths the fit gets worse. 

In conclusion we state that the thermal analysis of particle ratios allows for 
scaling of hadron masses. Such changes, however, result in modifications of 
the thermodynamical parameters for which the fits are optimal. In particular, 
lowering of all the masses leads to a smaller values of T chem and pf hem - This is 
a desired effect, since T chem ~170 MeV is large and may correspond to quark- 
gluon plasma rather than to a hadron gas. Our study of the modifications of 
the hadron widths shows that they have small impact on the ratios. 
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